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Here, a scenario is proposed, according to which a generic self-organized critical (SOC) system 
can be looked upon as a Witten-type topological field theory (W-TFT) with spontaneously broken 
Becchi-Rouet-Stora-Tyutin (BRST) symmetry. One of the conditions for the SOC is the slow driv- 
ing noise, which unambiguously suggests Stratonovich interpretation of the corresponding stochas- 
tic differential equation (SDE). This, in turn, necessitates the use of Parisi-Sourlas-Wu stochastic 
quantization procedure, which straightforwardly leads to a model with BRST-exact action, i.e., to 
a W-TFT. In the parameter space of the SDE, there must exist full-dimensional regions where the 
BRST-symmetry is spontaneously broken by instantons, which in the context of SOC are essentially 
avalanches. In these regions, the avalanche-type SOC dynamics is liberated from overwise a rightful 
dynamics-less W-TFT, and a Goldstone mode of Fadeev-Popov ghosts exists. Goldstinos represent 
modulii of instantons (avalanches) and being gapless are responsible for the critical avalanche distri- 
bution in the low-energy, long-wavelength limit. The above arguments are robust against moderate 
variations of the SDE's parameters and the criticality is "self-tuned". The proposition of this paper 
suggests that the machinery of W-TFTs may find its applications in many different areas of modern 
science studying various physical realizations of SOC. It also suggests that there may in principle 
exist a connection between some of SOC's and the concept of topological quantum computing. 

PACS numbers: 05.04.-a, 11.25.Tq, 12.40.Nn, 87.85.dm 



I. INTRODUCTION 

Self-organized criticality (SOC) [I'l is rightfully con- 
sidered by many to be a very fundamentalphenomenon. 
SOC found its applications in geophysics, [3 neuro- 
and evolutionary Q biology, cosmology and astro- 
physics, collective human (traffic flow, [8] stockmar- 
kets @) and animal behavior, cellular automation, 
[H in and many other areas of modern scientific re- 
search, [l^ [31 Previous investigations have firmly es- 
tablished several distinct conditions for and properties of 
SOC, [l^ which we discuss first. 

In Ref. SOC has been given another, yet more 

definitive name of .slowly driven, interaction dominated 
threshold systems. Expanding on this definition, SOC 
is observed in highly nonlinear systems, which possess 
large number of metastable states and which are driven 
by slow external stochastic noise. The noise pumps the 
energy into the system, which thus must also be capable 
of damping the excess energy, i.e, the system is dissipa- 
tive. What has been just said can be reformulated into 
the two following conditions for the SOCs that are thus 
expected to appear in 

(CI) stochastic nonlinear dissipative systems with 
large number of metastable states, when 

(C2) the external noise is "slower" than the internal 
processes (scale separation principle). 

SOCs possess a unique set of properties. The most 
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important of them are the following three: 

(PI) The time evolution is an infinite sequence of 
jumps between metastable states, i.e., of avalanches 
(P2) with algebraic correlations (criticality), 
(P3) which persist on moderate variations of the sys- 
tem's parameters (self-tuning). 

It is not known yet if there exist a mathematical con- 
struction that starts with conditions (CI) and (C2) and 
reproduces properties (PI)-(P3). In this paper, it is pro- 
posed that such construction does exist and is known 
under the name of Witten-type topological field theories 
(W-TFT). A generic SOC state must be identified as a 
W-TFT with spontaneously broken Becchi-Rouet-Stora- 
Tyutin (BRST) symmetry. 

The idea that a theory of SOCs must possess fermionic 
symmetries is not new. One of the known members of 
the SOC family is spin glasses. [15| The glass phase of 
spin-glasses, which is always critical, was identified [l6j 
as the supersymmetry (susy) broken state of an N=2 
susy theory, obtained from Parisi-Sourlas-Wu quantiza- 
tion of a Langevin stochastic differential equation (SDE). 
Moreover, the N=2 susy theory of self-organizing sys- 
tems has also been proposed, [l^ again based on the 
stochastic quantization of Langevin SDEs. At the same 
time, it is understood that a generic SOC system corre- 
sponds not to a Langevin SDE but rather to a general 
form SDE. Here, Parisi-Sourlas-Wu method is applied to 
general form SDEs, which is seemingly the only original 
part of the paper. For its most part, however, the pa- 
per is merely a collection of already established results, 
compiled into a physical picture, which we believe stands 
— behind all SOC phenomena. 

igmail.com 

We begin in Sec. |IT] with the discussion of a physi- 
cal difference between Stratonovich and Ito interpreta- 
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tions of SDEs. We show that the scale-separation prin- 
ciple of SOC unambiguously suggests Stratonovich in- 
terpretation. In Sec. Illli we demonstrate that the corre- 
sponding Parisi-Sourlas-Wu stochastic quantization leads 
to a (pseudo-Hermitian) model with a BRST-exact ac- 
tion, i.e., to a W-TFT as discussed in Sec. IIV Al In 
Sec lIVBl we discuss the localization principle and the 
need for a mechanism that breaks the BRST-symmetry 
and liberates SOC dynamics. In Sec. IIV C[ we show 
that this mechanism is instantons that in the context 
of SOC are essentially the avalanches. In Sec. lIVDi 
we speculate that BRST-symmetry breaking must oc- 
cur in full-dimensional regions of the SDE's parameter 
space. In Sec|Vl we switch to higher-dimensional SDEs 
and bring up the standard argument that the low-energy, 
long- wavelength part of the liberated SOC dynamics is 
the Goldstone Fadeev-Popov ghosts, which are gapless 
and thus are responsible for the critical distribution of 
avalanches (instantons). We argue that the "self-tuning" 
property is in the possibility to moderately vary the SDE 
parameters without hindering all the above reasonings. 
In Sec lVn we make a few remarks on the proposed sce- 
nario. We conclude in Sec lVIII 

II. SOC AS SDE 

A. Physical meaning of SDE 

Condition (CI) in the Introduction is essentially the 
statement that the natural starting point for the studies 
of SOC's are nonlinear dissipative SDEs. Consider an 
SDE for N stochastic variables: 

dt^\t)+A\t)^eit), (1) 

where i = 1...N, = A^{ip) is the vector field, which 
could be called the grift term, and is the stochastic 
noise. Mctastable states correspond to the critical points 
of the drift term. Consequently, for SDEs representing 
SOC systems there is a large number of critical points 
# {^a\A^{(fa) = 0} 1. This is yet another way of saying 
that the system under consideration is highly nonlinear. 
The drift term can also be given as 

= + (2) 

where the Langevin part can be defined via a potential, 
V, V ^ = 6^^V'j, V'j = dV/d(p^ (the summation over the 
repeated indices is assumed throughout the paper), and 
the non-potential (magnetic, Hamilton) part, A^, is such 
that A'^'^ = 5^^A]^. = -A^'^. For the arguments of this 
section it suffice to consider the Euclidian target man- 
ifold. In the following sections, we will generalize the 
discussion to Riemannian target spaces. 

That the system is dissipative suggests that V is non- 
zero. The potential part of the drift term is responsible 
for the tendency of the open SOC system to dissipate 
its energy into a reservoir and minimize the potential 



V , while the noise stochastically pumps the energy from 
(yet another) reservoir. The dynamics represented by 
Eq.(IT]) has the physical meaning of a stochastic energy 
fiow through a (highly non-linear) SOC system from an 
energy source to a drain. 

The stochasticity comes into the system only from 
the source. This assumes that we have already out- 
integrated the drain's degrees of freedom. A drain reser- 
voir with a memory can in fact introduce a tempo- 
ral non-locality into the SDE. In the lowest order ap- 
proximation this must have the following form: A' 
f_^M]{t - t')A>{ip{t'))dt' , where M could be called a 
memory kernel. The kernel has a certain characteristic 
time, Ad,M*(t — t') — > as \t — t'\/\d — > oo. In our in- 
terpretation of the scale separation principle, is (one 
of) the shortest time-scale(s) in the problem. Under this 
condition, one can always work on a scale much larger 
than Ad, on which the drain reservoir is memory- less and 
is rightfully represented as the Langevin part of Eq.([2]). 

Note also, that the "source" reservoir's noise, ^, may 
not necessarily be temporarily local (white). In partic- 
ular, our argument toward Stratanonvich interpretation 
of the SDE in the next subsection does not rely on the 
assumption of the "whiteness" of the source noise. 

For many (if not all) systems identified as SOC, the 
non-potential part of the drift term is also nonvanishing. 
Therefore, we are interested in cases when both V ^ Q 
and A' 7^ 0. Throughout the paper we refer to this situ- 
ation as to a general form SDE as opposed for example 
to the case of Langevin stochastic differential equation, 
A^ = 0, corresponding to Witten model, [ll] which is the 
(O-l-l)-dimensional beginning of all the W-TFTs. 

B. Stratonovich vs. Ito interpretations of SDE 

SDEs can be treated on field-theoretic grounds. At 
this, the quantization procedure depends on interpreta- 
tion of stochastic noise. There are two major choices: 
Stratonovich and Ito interpretations. They are related 
respectively to the stochastic quantization procedures of 
Parisi-Sourlas-Wu [li-[2l| and of Martin-Sigia-Rose. ^ 
The choice depends on physical conditions. In this sec- 
tion we show that the scale separation principle (condi- 
tion (C2) in the Introduction) unambiguously suggests 
the Stratonovich picture. 

Consider the discrete version of Eq. ^ : [2^ 

(^j+. - + ((1 - C)Aj+. + CAD = ?,V„ (3) 

where time now takes on discrete values separated by e, 
the noise on the rhs should be thought of as that acting 
during time interval [i, t + t], and C is the parameter that 
essentially represents the speed of noise in comparison 
with internal processes, e.g., to equilibration. Indeed, 
when the system is slower than the noise, it does not 
have enough time to adjust its variables to the quickly 
changing noise. Therefore, the infinitesimal change in <pt 
within each time step [t, t + e\ must be determined by the 
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earliest value of A% i.e., hy Al. Consequently, fast noises 
correspond to the Ito case with C = 1- 

In the opposite case of fast system, the Stratonovich 
choice of C = 1/2 is natural. The infinitesimal change in 
ipt is determined by " averaged" over the time- interval 
[t, t + e], or rather by its value in the middle of [t, t + e], 
i.e., by {Al^^^ + A\)I2. [Ill This essentially means that 
the system has enough time to adjust its variables to the 
noise, before the slow noise changes considerably. 

The only statistical partition function in the model is 
that of the noise: 



(4) 



where P is the noise distribution function. If we impose 
periodic boundary conditions on ip, however, the numbers 
of (^t's and ^t's will be the same. In this case the parti- 
tion function can be rewritten as a path-integral over ip^s 
instead of ^'s by the appropriate variable transformation: 



[d^]J(C)F(e(^)), 



(5) 



where ^{(p) is the Ihs of Eq.Q and the Jacobian of the 
variable transformation, J{() = \d£,t/dipt'\, is 



^(0 = 



Do 

No De 

D2, 







Nt 





DT^e 



with Dt 



Jt = + QA]^ and Nt = -t'H] -f C^?^-. 

Note that noise is typically assumed Gaussian and 
physics contained in Eq. ([4]) is trivial. On the other 
hand, Eq.® has appeared from Eq. ^ by the formal 
redefinition of the variables. Therefore, if the physics 
of Eq.® is not trivial this can only be blamed on the 
non-trivial topology {e.g., not one-to-one) of the hi ghly 
nonlinear (and nonlocal) map, (p : Ci^) ~^ ^^{i)- [IJ 
This is the first indication on the topological nature of 
the Parisi-Sourlas-Wu stochastic quantization procedure, 
the core of which as compared to Martin-Sigia-Rose pro- 
cedure is in "not forgetting" the Jacobian of the variable 
transformation. The physical justification for this appre- 
ciation of the Jacobian follows. 

On taking the continuous limit one obtains (up to a 
e-dependent constant) 



(6) 



For Ito interpretation, when C = 1, the first term in 
rhs of Eq.® is a constant. Furthermore, one can as- 
sume that due to the tendency of the system to minimize 
the potential, it spends most of its time in the region(s) 

where A\ = v'\ > 0. Therefore, J^o'^^i 



the second term is negligible, and one sets J(l) 1. 
Consequently, if the noise is fast one can neglect the Ja- 
cobian. 

The scale separation principle, however, suggests that 
for SOCs, the Stratonovich interpretation of stochastic- 
ity is more appropriate. Recall that in the sand-pile 
model, before one adds yet another grain to a random 
site, he has to time-propagate the system until it finds 
its new metastable state. Among the physical examples 
are earthquakes which are instantaneous when compared 
to the slow process of tension build-up in the earth crust 
due to plate tectonics, while the relatively long periods 
of (quasi-)equilibrium in the punctuation theory of bi- 
ological evolution are followed by sudden reconfigura- 
tions. Therefore, for the studies of SOC we must use 
Stratonovich interpretation of SDEs and keep the Jaco- 
bian with C = 1/2. 

One may expect, however, that for realistic situations 
C is not exactly 1/2. To see how C > 1/2 can change our 
story it is convenient to rewrite the Jacobian as: 



J{C) 



(7) 



where Jp — 2 sinh((l/2) J^^g AJJ can (and will be) 
represented as a path integral over additional fermions 
(ghosts) with periodic boundary conditions (see below). 
C > 1/2 explicitly breaks BRST-symmetry of the model. 
|26| The ghosts become energetically more costly due 
to the additional term {( — 1/2) Ai^ in the Lagrangian. 
This may lead to considerable changes in the low-energy 
dynamics even for small C — 1/2 ^ 1. Therefore, the 
Stratonovich interpretation of SDE is of cause an ap- 
proximation. It describes a hypothetical SOC system, 
for example a numerical model, for which the noise can 
be assumed infinitely slow. Having said that, we always 
assume C = 1/2 in the rest of the paper. 

So far we did not specify the noise weighting function, 
P. From previous discussion, a reasonable approximation 
for P of a slow noise is a temporarily non-local Gaussian 
with considerably large noise-noise correlation time, A^. 
At the same time, in the following sections we will use 
"white" noise as the driving force (As 0). This may 
seem a contradiction with the scale-separation principle. 
There is no contradiction though. One can straightfor- 
wardly generalize the developments in the following sec- 
tions to temporarily nonlocal noises, which, however, will 
turn out later to be an unnecessary complication. The 
point is that we can always rescale time t At so that 
As As /A 0. [13 From this perspective, the scale 
separation principle is solely in the appropriate choice of 
the stochastic quantization procedure. 
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III. PARISI-SOURLAS-WU QUANTIZATION 
OF SDE 

A. Path-integral approach 

Let us now proceed with the case of Gaussian white 
noise and with the Parisi-Sourlas-Wu quantization [2I] 
of Eq. (HI , apphcabihty of which is based on the appreci- 
ation of the importance of the Jacobian corresponding to 
the Stratonovich interpretation of SDE in the previous 
section. The partition function is: 

Z = l[d^]Jpe-floG^.^'^y\ (8) 

where — dt^p^ + and G*-' — {Gij)~^ is the noise- 
noise correlator, 

{e{t)e{t'))^G^'5{t~t'), (9) 

which for now is assumed to be independent of Lp (see 
below for the covariant generalization). One can always 
think of Gij as of the metric of the target manifold as long 
as is arbitrary. In other words, there is no Riemannian 
structure on the target manifold yet, and © could as well 
play the role of the metric. It will be seen later that it 
actually does play this role. 

The Jacobian can be represented as the path inte- 
gral over the fermionic Fadeev-Popov ghosts according 
to detAf — /[dx][dx] exp(xiMjx-' )- Furthermore, one 
can employ the Legandre multiplier, B^, which is the dy- 
namical conjugate to Lp. The partition function now is: 




where $ represents all the fields, and the action S = 
p L is defined by the Lagrangian 

L = iB.{lC - iG''B,/2) - Udtx' + A^,x')- (H) 

The periodic boundary conditions are imposed on all the 
fields. In case of ghosts, these conditions are needed if 
the fermion determinant is to represent the Jacobian. 

The model enjoys the global nilpotent fermionic 
BRST-symmetry (Q-symmetry) [2^ given by the in- 
finitesimal operator 

2 = E- / x\t)S/Sip\t)+tB,{t)d/dMt), (12) 

so that QS = {Q, S} — 0. Importantly, the action is also 
Q-exact: 

5 = {Q,*}, (13a) 

where 

* = r UK'-iG^'B,/2), (13b) 
Jt=o 



is known as the gauge fermion. Q-exact action is a unique 
feature of W-TFTs. It looks like the whole action is 
nothing else but the BRST gauge-fixing term. In Ref. [s^] , 
this was identified as the "quantizing zero" situation. 

From the field-theoretic point of view, a more appeal- 
ing derivation of is based on the philosophy that the 
SDE itself can be looked upon as the gauge choice for the 
two "independent" stochastic variables - if and ^: 

Z = J [d^][de][dx][dx]e-^-oG..«'«V2+...^ (14) 

where dots denote the gauge-fixing Q-exact term, 
{2Jt=o^»(^' - r)} (here Q has the form (HH with 
iB, ^ G,,^^). By noticing that G„e^V2 = {Q, M'J^} 
is also Q-exact and by formal redefinition, ^' ^ iG'^^Bj, 
one recovers (fT3|) . 

Moreover, the above Parisi-Sourlas-Wu quantization 
can be generalized to the cases when the noise-noise cor- 
relator (the metric of the target manifold) is dependent 
of G*-' —7- G'^^{(p). The Batalin-Vilkovisky procedure 
must be used to come up with the same Q-exact action 
([T^ . [lilPl At this, however, the Q operator ^ will 
acquire a more intricate form, which accounts for the 
curvature of the target manifold. The detailed explana- 
tion of the corresponding quantization procedure and the 
appropriate form of the Q-operator can be found in Ap- 
pendix A. 2 of Ref. [30j . where the only adjustment to our 
case needed to be done is G^^Vj — ?• A^ . 

After the covariant generalization, Eq. (jl3p represents 
a very general (0-1-1) SDE, which is quantized stochas- 
tically by the Parisi-Sourlas-Wu method in accordance 
with the Stratonovich interpretation of the slow noise. 

B. Schrodinger picture 

Let us turn now to the Schrodinger picture. The sta- 
tistical (Euclidian) partition function ([TU)) can be given 
through the Hamiltonian function as : 

Z = J [d$]e^*=o(*"^'^"^"+"x'«'x'-^), (15) 

where tt's are the canonical momenta, which on passing to 
the Schrodinger picture w^i ~> —id^i and vr^i — > —id^i. 
From Eqs.([T5|) and (|lip we identify n^i = —Bi and vr^i — 
iXi, so that in flat coordinates: 

H = -A/2-[A^9^>] + /2 + AJ^.[a^.,x^■]-/2, (16) 

where A = d^^iG^^d^j is the Laplacian. The choice of 
the operator ordering in the second term of Eq. (jl6l) is 
standard. The operator ordering in the last term has 
direct connection to the Stratonovich interpretation of 
the noise of the SDE. Have we been considering "faster" 
noises, we would have to use (A^jd^iX'' ~ (1 ~ O^^'jX^d^i 
instead of the last term in Eq. ([T5)) , that is we would have 
to add (C - 1/2) A?, to the Hamiltonian. (c.f., Eq.©) 
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Noether charge associated with the Q-symmetry is: 

Q = -iB,x' ^X'd^., (17) 

which is nothing else but the exterior derivative on the 
target manifold. This is so far the second explicit reveal- 
ing of the topological nature of the model. As it should, 
the charge is nilpotent, = 0, and commutative with 
H. 

The Hamiltonian can also be given as: 



H = {Q,Q}/2, 



(18) 



where 



X^{iG''BJ + 2A') ^ d^. i-G'^d^, - 2A'). (19) 



In non-flat coordinates, the Hamiltonian is again an an- 
ticommutator (|18p. while explicitly: 



H 



-A/2 - Ca^ 



where the Laplacian is given by the Weitzenbock formula 
that includes four-fermion coupling through the Riem- 
mann curvature tensor, Cai is the Lie derivative along 
A*. Noether charge is the same, while Q = *Q * — 2i(A*) 
with ★ and i(A*) denoting the Hedge operation and inte- 
rior multiplication by A^. 

Hamiltonian (|20p has a very clear physical meaning. 
The first term is the quantum mechanical "smearing" 
(dispersion) of wave-functions, which are p-forms from 
the (complexified) exterior algebra of the target mani- 
fold. The intensity of the dispersion is determined quan- 
titatively by the "magnitude" of noise- noise correlator, 
I |G*^ 1 1 ^ temp, which thus has the meaning of noise tem- 
perature. This is the essence of the stochastic quantiza- 
tion, the stochasticity takes the form of the quantum me- 
chanical fiuctuations. In the classical, low-temperature 
limit, [slf one is left with a non-dispersive classical flow 
of p-forms along the drift term, A*. 

Hamiltonian (|20l) is not Hermitian with res pec t to the 
conventional metric in the Hilbert space: [sol |32| 



(*a*)A/3, 



(21) 



where A is wedge product of p-forms and the integration 
is over the target manifold. All entries in the Hamilto- 
nian, however, are real. Therefore, the Hamiltonian can 
be looked upon as an infinite-dimensional real matrix. 
The spectrum of such Hamiltonian consists of real ener- 
gies or pairs of complex-conjugate energies. This means 
that the Hamiltonian is pseudo-Hermitian [s^ and there 
must exist such hermitian, invertible rj that 



= rjHrj 



(22) 



Model ()13p is a pseudo-Hermitian quantum mechanics 
[33„ ^4] with T] being the metric of the Hilbert space 

{{c^m = (^hv), ((-^i = (0k/, 1^)) = iv-), (23) 



preserved by the Schrodinger evolution: 

^dtMyj}} = MH - T]-^H^Tim = 0. (24) 

For Witten model, rjw — . In fact, Witten model is 
a quasi-Hermitian quantum mechanics. The eigenvalues 
of its Hamitonian are all real, since rj^^^ Hrj^^^'^ (such 
transformation brings Eq.(20) to its conventional form 
appearing in the literature) is Hermitian. For a general 
form SDE, 77 is complicated and highly non-local, and the 
energies are complex. 

The anticommutator form of the Hamiltonian (IT51) is 
that of the N=2 pseudo-supersymmetric (p-susy) quan- 
tum mechanics, 35] with the following operator algebra: 



i^o = {Qo,Q&}/2,Q^-[i/o,Qo]=0, 



(25a) 
(25b) 



(20) 

and Q\) = rj ^Q^rj is the pseudo-Hermitian conjugate to 



Qo- In order for the p-susy to be unbroken, the opera- 
tor algebra (|25|) must be complemented with the ground 
state(s) such that 



go|o)) = Q«|o)) = o. 



(26) 



N=2 p-susy is twice larger than the Q-symmetry. It is 
a combination of Q-symmetry and yet another fermionic 
pseudo-anti-Q symmetry related to Q by the pseudo- 
time-reversal conjugation (777~-conjugation). 

In general case, however, model (fT3|) does not possess 
the pseudo-anti-Q symmetry. Indeed, Q, which is sup- 
posed to play the role of Qq, is neither nilpotent nor does 
it commute with H. With Q for Qq, the operator algebra 
Eas. (j25bp fails. For example, in flat coordinates: 



(g)2 = 2F^G'=^a^,a^. ^0, 



(27) 



where = {AJj - A^^)/2 is the "field tensor" of the 
non-potential (magnetic) part of the drift term. 

If for Fj ^ the N=2 p-susy is present, it means that 



(g)2,[i/,g]|phys)) = 0, 



(28) 



are satisfied for all the (physical) states of the model. For 
these situations, one can find such ij that Q = rj~^Q'^-q. 
These situation may appear, for example, if A^ is a Killing 
or a symplectic vector field, [s^, [13] There is no reason, 
however, to believe that this is true in general. Below 
we will argue that in the SDE's parameter space there 
must exist region(s) with explicitly broken N=2 p-susy 
(together with the spontaneously broken Q-symmetry). 

In the rest of the paper we use the path-integral repre- 
sentation of the pseudo-Hermitian quantum mechanics. 
This is advantageous for the reason that finding the ex- 
act form of the metric, 77, is a complicated task and one 
should avoid rj whenever possible. Fortunately, the met- 
ric is automatically incorporated into the path-i nteg rals 
and one must not worry about its explicit form. 
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IV. SDE AS W-TFT 
A. Conditions for W-TFT 

To identify a theory as a W-TFT, one needs: [s^l 
(i) a Z2 -graded nilpotent fermionic Q, 
(m) a Q-exact action, 

{Hi) a Q-invariant path-integral measure, and 
(iv) a Q-closed ground state(s), [s^ which ensures that 
the Q-symmetry is not spontaneously broken. 

These conditions suffice to establish a unique set of 
properties of W-TFTs. [13] In particular, one can intro- 
duce a "metric" for time: dt e{t)dt, and in case of 
higher-dimensional theories the metric for spatial dimen- 
sions, g (our case is lacking Lorcntz invariance so that e 
and g should not be combined into a space-time metric) . 
In result, the gauge fermion acquires explicit dependence 
on e and 5: ^ ^ '^{e,g). The topological nature of W- 
TFTs is seen through the (e, g)-independence of (now we 
consider T— >■ 00 limit corresponding to the field-theoretic 
interpretation of the path-integral): 



mrAB. 



.|0)) = J [d$]e-{2,*}^5. 



(29) 



where A,B... are Q-closed, i.e., {Q,A},{Q,B}... = 0, 
and T denotes chronological ordering. Indeed, the func- 
tional variation with respect to, e.g., e{t): 

mTAB...\0)) = {{0\T{Q,iemAB...\0)) = 0(30) 



6e{t) 



where 7e(i) = —S'^/5e{t). The last equality can be 
proven by partial integration (Q is a differentiation op- 
erator) and holds only if the ground state(s) obey (in the 
Schodinger picture) 



Q|o)) =.g«|o)) =.0, 



(31) 



with Q'^ = rj ^Q^rj. We see that this condition (and/or 
condition (iv) above) is actually that for the unbroken 
N=2 p-susy (HH). 

We can not speak with certainty that the breakdown 
of N=2 p-susy always tailors the spontaneous breaking of 
Q-symmetry. There may exist situations when (5|0)) =0, 
while Q'lO)) 7^ 0. Such situations are at least exotic. We 
omit these exotic situations in this paper and think of 
the two symmetry breakdowns as of equivalent. 

In fact, the proof of that the Q-symmetry is unbroken 
is a subtle and comphcated part of W-TFTs. [H E^l 
Therefore, we find ourselves now in a fortunate position 
that the Q-symmetry breaking is actually what we are 
looking for. To demonstrate the breakdown of the Q- 
symmetry of (|13l) it suffice to establish, for example, the 
T-dependence of the partition function. Indeed, the par- 
tition function is one of the topological invariants ((29)) 
with AB... being a unity operator that is obviously Q- 
closed. This suggests that Z of a true W-TFT is time- 
deformation invariant and thus T- independent. On the 



other hand, model ([13]) definitely satisfies conditions (z)- 
{iii) above so that the spontaneous Q-symmetry breaking 
is the only possible source for the T-dependence of Z of 
otherwise rightful W-TFT with T-independent Z. 



B. Localization principle 

Let us turn now to the (lowest order) one-loop study of 
partition function. Consider the case of low temperatures 
(I |G'-' II ^ 1) so that A^ is very "pronounced" . The major 
contribution into the partition function comes from the 
lowest energy states that cluster around critical points of 
A (the metastable states of the SOC). This can be seen 
by out-integrating the Legandre multiplier, B: 



Z = / [d*] 



- f^^g K'GijK' /2+formions 



(32) 



The contribution is that from the fluctuations around the 
stationary paths, (p*i(i) = so that: 



_2:(i-ioop) = ^ 



(33) 



Za can be calculated in the locally Euclidian coordinates, 
5if^ , in which the Gaussian part of the action is: 



^(2 



(5^\'D^V+r5^^l2 + x'fflx') , (34) 



where 2?+ — dt+a, P_ = —dt+a^ with a^^ being A'^ ((/?„) 
in the (5<y9-coordinatcs. The Gaussian integration leads to 
the bosonic and fermionic determinants: 



(35) 



Introducing temporal Furrier components, e.g., Stp'' = 
En=-oo S^'in)e'^-\uj^ = 2^7i/T, one gets 



X a / aa 



T|l/2 



n 



j:J-^i\ujl + a^\\Lol + {aT)2\y/2- 

To simplify further this expression one notices that since 
a is real, its eigenvalues, ai,i — 1...N, are either real 
or come in complex conjugate pairs, and that the set of 
eigenvalues of a'^ = is a*,i = 1...N. Therefore, the 
sets of eigenvalues and consequently the determinants of 
-I- and uj^ + (a^)^ are the same. Hence, Ea — I 
and: 

^a = (-l)^", 

where Aa is the number of real negative e„'s. (—1)^° = 
sign|A*((/3Q)| = inda is known as the index of the critical 
point, and 



^(i-ioop)^y ind„, 



(36) 
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is the Euler characteristic of the target manifold accord- 
ing to the Poincare-Hoft theorem. Hence, on the one-loop 
level, the fluctuations leave the Q-symmetry intact as is 
seen from the T-independence of Eq. ([55]) . 

The reason for bringing the reader's attention to the 
one-loop calculation that does not break the Q-symmetry 
is twofold. Firstly, Eq. ((36)) is yet another explicit demon- 
stration of the topological nature of the model. Sec- 
ondary, it reveals the tendency of ghosts to compen- 
sate for bosonic fluctuations. Up to one-loop, the ghosts 
completely cancel fluctuations in the bosonic fields which 
means that so far the model has no dynamics. This com- 
pensation is known in the literature as the localization 
principle. As the name suggests, the path integral is "lo- 
calized" to classical solutions of the SDE. Stationary so- 
lutions at critical points is one class of classical solutions. 
The other class is instantons (see next subsection). 

From the physical point of view, the meaning of the lo- 
calization principle is that the slow noise does not provide 
the system with "frequencies". In result, the dissipative 
system does not fluctuate. 

For Witten model, the one-loop approximation is 
known to be exact from the perturbative point of view. 
For a general form SDE, however, higher-order fluctua- 
tional corrections may in principle break the Q-symmetry 
(anomalously in that sense that perturbatively provided 
corrections possess lower symmetry than the original ac- 
tion). [4l| However, even if the fluctuations leave the 
symmetry intact, it can still be broken by instantons. 
[so, .32. ,4Q. i43i| In fact, instantons is the primarily source 
of the Q-symmetry breakdown as we discuss in the next 
section. In particular, instantons break the Q-symmetry 
even of Witten model itself, when the Euler characteris- 
tic of the target manifold is zero, (s^l It is the instanton- 
induced Q-symmetry breakdown which we identify as an 
SOC and which we are interested in. 



C. Meaning of instantons and their role in 
Q-symmetry breaking 

The SOC dynamics can be described as follows. The 
system spends most of its time in its metastable states, 
while the time evolution is an infinite sequence of sud- 
den jumps (avalanches) between the metastable states. 
In our terms, the metastable states are nothing else but 
the perturbative ground states (PGS), \a)), around crit- 
ical points of A^. The "jumps" between different PGS's 
are the processes of the quantum mechanical tunneling, 
which are called instantons. Therefore, in the application 
to the SOC systems, instantons are nothing else but the 
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mathematical term for the avalanches. 

Instantons in a sense are the opposite of fluctuations. 
The processes of quantum mechanical tunneling are ex- 
ponentially weak as compared with fluctuations. That 
instantons are of low-energy means that they do not 
happen often. Fluctuations are localized around criti- 
cal points unlike instantons which connect them. Hence, 




FIG. 1: (Color online) (a) Schematics representing an SDE's 
drift vector field, A^, and its critical points (circles) with 
their indices given explicitly. Thick double- arrowed curves 
represent instantons (classical solutions of SDE connecting 
two critical points). They are complemented by the axil- 
lary thin vector-lines which connect critical points with in- 
dices differing by 2. Thick "reversed" double-arrowed curves 
(given for only one closed path) represent antiinstantons, 
which as is explained in the text have different matrix el- 
ements from those of the corresponding instantons due to 
the non-potential (magnetic) part of A^. This imbalances 
instantons and antiinstantons on time-reversed loops, e.g., 
the inner and outer loops shown. Under certain conditions, 
this imbalance must lead to instanton-induced spontaneous 
BRST-symmetry breaking. The argument is similar to that 
of Anderson localization of a particle on a random lattice. In 
the presence of a magnetic field, the constructive interference 
of time-reversed hopping loops, which is responsible for the 
localization, fails, (b) In the semi-classical approximation, 
each instanton is a Gaussian path-integral around classical 
solution of the SDE (thick center curve) connecting two crit- 
ical points {ipa and if^) ai t — ±oo. This solution has the 
so-called modulus - its center. The variation in the modulus, 
e.g., ipci{t) — >■ (pci{t ± St), does not change the action of the 
instanton. This leads to the existence of a zero fiuctuational 
mode, 5ipo{t) = dtfciit). The corresponding zero ghost mode 
decouples from the action and the path-integral is zero unless 
this zero mode is matched. The matrix element of Q-operator 
is non-zero because it also has an additional ghost operator 
matching the zero ghost mode. By similar arguments, instan- 
tons can only connect critical points with indices differing by 
1 (in case of antiinstantons, —1) as is seen from figure (a), 
(c) The qualitative difference between excitation spectra of 
higher-dimensional theories with unbroken (left) and sponta- 
neously broken (right) Q-symmetry. The double curve on the 
left represents dispersions of supersymmetric partners. The 
cone on the right represents goldstinos. The negative energy 
states (shaded) are occupied and form a (gap-less) Dirac sea. 
The Dirac-sea ground state represents a quantum liquid of 
solitons (textures, patterns) as explained in the text. 



instatons are of ultimate nonlinear character. In linear 
situations they simply do not exist. This is yet another 
sense in which instantons are qualitatively different from 
fluctuations. Furthermore, as it will be clear later from 
the discussion of higher-dimensional models, avalanches 
are the only low-energy, long- wavelength dynamics in the 
system. This is especially appealing for the physical pic- 
ture of SOCs. 

The above picture of SOC makes sense only when the 
critical points possess distinct PGS's. In other words. 
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when the PGS's do not overlap much. The overlap can 
be provided by a sufficient "smearing" of the PGS's due 
to large kinetic energy. In the model under consideration, 
the kinetic energy in (|20p is determined by the noise tem- 
perature (the magnitude of the metric G"-' ). Therefore, 
for an SOC to occur (at least in its recognizable form) 
the noise must not only be slow but also weak. 

The low-energy instanton dynamics can be defined by 
the projection onto the reduced Hilbert space spanned 
by |a))'s: [Mill 

H«;3 = {Q,Q}a^/2, 

= {{a\Q\f3)),Qa.0^{{a\Q\f3)), 

where Q and Q are matrix elements of instantons and 
antiinstantons. In the path-integral language, 

\a)) = j [d^e- ^ , (37a) 

where the functional integration is over the paths starting 
at t = at the argument of |a)) and ending at i = +oo 
at iy9^. The (pseudo-) bra is accordingly: 




(37b) 



where the integration is over the paths starting a,i t — 
— oo at (/?^ and ending at t = at the argument of {{a\. 

From Eas.(l571). the instanton and antiinstanton matrix 
elements can be given as: 

Qcp = y"[d$]Q(0)e-^*=--^, (38a) 

Qa/3 = y"[d$](5(0)e-^*i-~^, (38b) 

where the integration is along the paths connecting the 
two critical points, (^*(— oo) = ip"^^ and (p^{+oo) = (/J^. 

In the one-loop approximation, the matrix elements of 
instantons can be found by Gaussian integration around 
a classical solution connecting the two critical points, 

dt^'^,{t) + A^{^,m = 0, = = 

As in Sec. (jIVB[) . the Gaussian integration can be done 
in the locally Euclidian coordinates for Lp, so that: 

Qa/^ = y"[d$]Q(0)e-^*'\ (39) 

with 5''^-' from Eg. (1341) . This time, however, oper- 
ators 'D± have explicit time-dependence, e.g., P+ = 
dt + a((/9*[(t)). If the classical solution connects critical 
points with indices differing by one, inda = ind/? -I- 1, 
operator "Dj^ has one zero mode. Indeed, due to the 
time-translation invariance, the classical solution can be 
shifted in time (p* j {t) j [t + 5t) and still remain the 
classical solution. Differentiation of SDE with respect 
to 5t results in V+dtp^yit) = 0. This mode corresponds 
to the infinitesimal shift of the center of the instanton, 
which is called instanton modulus. Furthermore, there is 



the corresponding ghost mode which decouples from the 
fermionic action. Therefore, the path-integral contains 
an integration over an unmatched Grassmann number. 
Such path integral must be zero unless the matrix ele- 
ment is that of an operator of ghost number 1. Q is such 
an operator and its matrix element is non-zero. 

Consider now the classical solution which connects crit- 
ical points with indices differing by more than one, e.g., 
inda — ind/? + 2, represented by thin lines in Fig.(IT])a. 
Such solution can be though of as an instanton followed 
by yet another instanton. Each instanton comes with 
its own modulus and with the corresponding zero ghost 
mode. The path-integral has two unmatched Grassmann 
numbers. The matrix element of the ghost number 1 
operator Q is zero. Hence, 

^ 0, inda = ind/3 + 1. (40) 

Similar reasonings apply to antiinstantons, which cor- 
respond to the backward propagation in time. Gaus- 
sian integration must be performed around a classi- 
cal solution of the "time- reversed" differential equation, 
dtvl,{t) - A\p,m = 0, (^^i(-oo) = (^^i(+oo) = v'^. 
The antiinstanton analogue of ([M)) is 

= g-2AV,„-2/A-G.,<i^^,, y"[d$]Q(0)e^^''\(41) 

where ^Vfia = V'(<P/3) — '^Wa) > and A is the non- 
potential part of A (see Eq.((2])). Now, it is D_ which has 
a zero mode and the corresponding x must be matched, 
which Q{0) actually does so that: 

7^ 0, inda + 1 = ind/3. (42) 

As compared to Ea. ((39|) . Eq. (|4T|) has two additional ele- 
ments. The first one is AVa„, which has the meaning of 
the potential difference at two critical points. It provides 
an antiinstanton with the exponentially weak tunneling 
factor. The factor can be absorbed into the pseudo- 
Hermitian metric of the reduced Hilbert space by the ap- 
propriate rescaling of PGS's: |a)) e^('^°'|a)), ((a| 
((ale"^'"^"^. In Witten model, this will result in that 
Qa0 — Q/3a- The other part, / A''Gijd(p-'^^, has the mean- 
ing of the (imaginary) phase factor acquired by a particle 
moving in a vector potential, A^, of a magnetic field. Due 
to this part, the matrix elements of instantons are differ- 
ent from those of the corresponding antiinstantons even if 
we try to rescale |a))'s. This imbalances instantons and 
antiinstantons. The hopping (instanton-antiinstanton) 
evolution along A^ becomes preferable. 

Importantly, the integrant in Eq. ([M|) is Q-exact: 
g(0) = -iB,{0)xHO) = -{Q,UO)x'm. a non-zero 
average of a Q-exact operator, is a direct indi- 

cation on the Q-symmetry breakdown. Therefore, the 
Q-symmetry is prone to be broken by instantons. (43j 
Thus, instantons (antiinstantons) is the primarily candi- 
date for the Q-symmetry breaking. In fact, this effect is 
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so pronounced that even balanced instantons and antiin- 
stantons of Witten model break Q-symmetry if the Euler 
characteristic of the target manifold is zero. [3C| 

The identification of as a vector potential of a mag- 
netic field brings about the analogy with the problem of 
Anderson localization (see FigH]). Originally, 4^ An- 
derson considered a quantum particle which lives on a 
random lattice (the "lattice" of the PGS in our case). 
The particle hops between the lattice sites (instantons 
and antiinstantons in our case). Due to the constructive 
interference between time-reversed paths, the probabil- 
ity to stay at the same site is always greater than that of 
traveling and if the conditions are right this leads to the 
localization (preservation of Q-symmetry in our case). 
Once, however, one introduces the magnetic field into the 
system, the time-revered paths get imbalanced and the 
constructive interference argument fails thus leading to 
the delocalization (Q-symmetry breakdown in our case). 

The conclusion of this subsection is as follows. If in the 
low-temperature limit the Q-symmetry is broken, it is 
most likely due to the instanton-antiinstanton imbalance 
induced by the non-potential part of the drift term. The 
low-energy part of the liberated dynamics is of instan- 
ton (avalanche) type and from the mathematical point of 
view is described by the low-energy ghost modes, which 
in turn represent modulii of instantons (avalanches). 

From the quantum-mechanical treatment, however, it 
is not clear why the distribution of avalanches (instan- 
tons) must be a power-law. To see this one must turn 
to higher-dimensional theories as we do in Sec|Vl Be- 
fore that, however, a few words on the conditions for the 
Q-symmetry breakdown are in order. 



D. BRST-symmetry breaking 

An important question is when the spontaneous break- 
down of Q-symmetry occurs. We do not have a satis- 
factory answer to this question. We believe, however, 
that the answer may have a lot to do with the interest- 
ing observation of Ref. [i^. There, it was suggested 
and demonstrated with a few examples that in VT- 
symmetric pseudo-Hermitian quantum mechanics the ex- 
plicit breakdown of N=2 p-susy is always accompanied by 
a spontaneous breakdown of PT-symmetry. The break- 
down(s) occurs when the parameters of a model reach 
some critical values and the pairs of complex-conjugate 
energies appear. PT-symmetric models [3J| is a subclass 
of pseudo-Hermitian jyT-symmetric quantum mechanics. 
|33| Therefore, it is reasonable to expect that the proposi- 
tion of Ref. [4^ can be generalized to other TyT-symmetric 
models such as the one under consideration. 

That the ry7"-symmetry breaking is needed in our case 
can be seen again from the analogy with the Anderson 
localization problem. In order to delocalize an Anderson 
particle one needs to break the T-symmetry (to imbal- 
ance 7~-reversed paths) by the introduction of a mag- 
netic field, whereas in our case one needs to break the 



77T-symmetry by the introduction of the magnetic (non- 
potential) part of the drift term. 

In our case, there is also a seeming indication on the 
necessity of the appearance of complex energies for the 
existence of the liberated dynamics. The point is that 
the Parisi-Sourlas-Wu stochastic quantization leads to a 
" statistical" (Euclidian) partition function. The operator 
of the quantum mechanical time-evolution has the form 
of the propagation in imaginary time: 



f ^ r. 



e-*^in))((n| 



(43) 



where n numerates the levels and En are their ener- 
gies. Imaginary time usually has the meaning of tem- 
perature. On the other hand, we know that the time in 
the evolution operator is the original time of the SDE 
and not a temperature of any kind (temperature in our 
case is the metric of the target manifold). Therefore, the 
time-evolution does not make a conventional quantum- 
mechanical sense. Presumably, this can be attributed 
to the absence of the propagating modes and/or to the 
unbroken Q-symmetry. 

When a pseudo-Hermitian Hamiltonian, however, pos- 
sesses complex energies. En ^ En + lEn, the evolution 
operator becomes: 



^ — 



(44) 



which is now more of the taste of a conventional quantum 
mechanical evolution. One can hypothesize that this is 
the result of the liberated dynamics and/or of the spon- 
taneous breakdown of the Q-symmetry. 

The emerging picture is as follows. If we take a Wit- 
ten model with unbroken Q-symmetry and start chang- 
ing the SDE's parameters, we will eventually reach some 
critical values at which Q- and yyT-symmetries will spon- 
taneously breakdown. This will be signified by the ap- 
pearance of complex energies of the Hamiltonian. 

Disregard of whether the speculations in this subsec- 
tion are correct or not, it is natural to expect that SOC 
dynamics is liberated from otherwise a rightful W-TFT 
in some low-temperature full- dimensional regions of the 
SDE's parameter space. At least for N=2 susy models of 
spin-glasses this statement seems to be correct, and 
there is no reason for the situation to be dramatically 
different in the more general case of N=2 p-susy. 



V. HIGHER-DIMENSIONS AND GOLDSTONE 
CRITICALITY 

The fact that the Q-symmetry is spontaneously bro- 
ken suffice to establish the criticality of the avalanche 
dynamics. All what is said for (0+l)-theories can be 
generalized to higher dimensions. [46j The generaliza- 
tion can be viewed as a limit of an infinite-dimensional 
target manifold. Literally, the index i is split into the 
spatial coordinates, x, and the coordinates of the target 
manifold: (yS* f^{x). 
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The higher-dimensional counterpart of Eq.([T]) is: 

dt^'{xt)+A'{^)^^'{xt), (45) 

where A*(v?) is some functional of (p^s and the Gaussian 
noise correlates on the metric of the target manifold: 

{Cixt)Cix't')} = g'^iip)6{t - t')5'^{x ~ x'). (46) 

After the Parisi-Sourlas-Wu quantization, the action is 
Q-exact and defined by the gauge fermion (|13bp with the 
additional integration over the spatial coordinates, x. 

In situations when the Q-symmetry is spontaneously 
broken there must exist a local operator p{xt) such that: 

{{{Q,p{xm^O. (47) 

By the standard arg ument, there is a gapless Goldstone 
ghost mode. [13, Consider the following average: 

{{p{xt))) = I mp{xt)e-''. (48) 

Under the path-integral we can make the space-time 
dependent infinitesimal transformation of the fields 
5,^{x't') = e{x't') {Q, ^{x't')}. The average is not to be 
changed since the transformation is merely the change in 
the integration variables: Se{{p{xt))) = or 

eixt){{{Q,pixt)})) = {{pixt)S,S)). (49) 

By the Q-exactness of the action: 

S,S^ f {{d^e)jnix't'), (50) 

Jx't' 

where index p combines space-time coordinates and J** is 
the Noether current associated with Q. e can be chosen 
constant within some arbitrarily large space-time volume, 
il, such that (xt) e il, and e = outside the volume. 
From Eqs.([50]), dlH), and (gT]) we get the integral form of 
the failed Ward-Takahashi identity: 

/ {{p{xt)J''{x't')))n{x't')^^ const ^0, (51) 

Jdn(x't') 

where the integration is over the boundary of il and 
is the unit vector normal to the boundary. As long as 
f2 can be chosen arbitrarily large, the above equation 
shows that {{p{xt)J^{x't'))) falls off algebraically as a 
function of cc' — a; and t — t' . This happened due to the 
existence of the gapless Goldstone ghosts (goldstinos), 
which are "struggling" to restore the broken symmetry 
of the vacuum. 

The spectrum of excitations must look like the one 
given in Fig[TJ;. Goldstinos must form a gapless Dirac 
sea with negative energy states being occupied. Fur- 
thermore, the would-be bosonic superpartners of goldsti- 
nos (upper curve of FigHJ;) must also be liberated from 
the fermionic symmetry and in many cases of interest 
they must be gapped. Thus, goldstinos (representing 



avalanches, see below) is the only dynamics the system 
has in the low-energy, long- wavelength limit. 

In fact, p{xt) in Ea. (|T7)) is the spatial density of ghosts 
and due to the existence of the Dirac sea, {{p{xt))) ^ 
0. In Ref . [l7l| . {{p{xt))) was identified with the entropy 
density. In this manner, the spontaneous Q-symmetry 
breaking and the appearance of the corresponding Dirac 
sea is related to the concept of " entropy production" . 

As in case of supersymmetric quantum mechanics, (32j 
the ghosts of PGS's (which in higher-dimensional cases 
could also be called perturbative vacua) represent unsta- 
ble directions of in the functional space of 'p{x). To 
find the ghosts of a particular perturbative vacuum, one 
must diagonalize the linear operator {5A^I5'^{x')') at 
the corresponding static field configuration, '/^^(a;), such 
that A^i^Lpa) = 0. Those modes which have negative real 
parts of their eigenvalues are occupied by the ghosts. 

In cases of interest, the drift term functional A^{ip) 
does not explicitly depend on the spatial position x. 
Therefore, the fact that a perturbative vacuum has many 
unstable directions (ghosts) most certainly implies that 
ipaix) is spatially inhomogeneous. Such vacua could be 
looked upon as textures, patterns, etc. They can also 
be called configurations of solitons such as domain walls, 
vortices, etc. The actual "Dirac sea" vacuum is a quan- 
tum mechanical superposition of perturbative vacua cor- 
responding to various solitonic configurations. In other 
words, the ground state (s) of an SOC is a quantum liquid 
of solitons (textures, patterns). 

As compared to the case of quantum mechanics dis- 
cussed previously, higher dimensional theories must have 
an important new element. Critical points of the func- 
tional are not isolated. [48|] Solitons could be thought 
of as instantons in space (not time) . Therefore, they also 
have modulii at least due to the symmetries of the space. 
For example, if the model is invariant with respect to 
translations in space, field configurations obtained from 
a given ipl^(x), A'^{(pa) = by all the spatial translations, 
faxi^) = + also critical field configura- 

tions for A^, A'^{^pax) — 0, and thus are also candidates 
for a perturbative vacua. In this case, Bott-Morse the- 
ory applies, js^ which states that the perturbative vacua 
are from the cohomology of the modulii space of solitonic 
configurations. For the soliton modulii of spatial transla- 
tions, this corresponds to the zero-momentum vacuum. 

As we discussed before, the low-energy, long- 
wavelength part of the liberated modes, which as we 
saw are the goldstinos, must represent modulii of in- 
stantons (avalanches) connecting different perturbative 
vacua. Therefore, the algebraic correlator of the gapless 
ghosts assumes critical distribution of avalanches. We 
believe this is the essence of the criticality of SOCs. 

From discussion in Sec lIVD] we know that in the 
parameter-space of the SDE, the model has Q-symmetry 
spontaneously broken in regions of the same dimension- 
ality as the parameter space itself. Therefore, for any 
SOC we have a freedom to moderately vary the SDE pa- 
rameters leaving the system at the Goldstone criticality. 
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We believe this is the essence of the self-tuning property 
of SOCs. 



VI. DISCUSSION 

In this section we would like to make a few remarks 
which seem interesting: 

• Previous studies of SOC used Ito interpretation of 
SDEs. This approach led to the conclusion that SOCs 
are members of the family of non-equilibrium phase tran- 
sitions such as directedpercolation [4§| captured by the 
Reggeon field theory. p(f\ 

Such approach certainly fails to explain the self-tuning 
property of SOCs. The point is that on the phase- 
diagram, i.e., the SDE's parameter space, a conventional 
critical state separates phases of different qualities and 
thus occupies manifolds of lower dimensionality than the 
phase-space itself. In other words, there is always at 
least one direction in the phase space, which leads the 
system off its criticality. The picture proposed in this 
paper seems to resolve this issue by explaining the crit- 
icality of SOC by the Goldstone theorem. Note, that 
the Goldstone explanation of the "self-tuning" property 
of s pin- glasses is known for both the N=2 susy approach 
[m Il6| and replica trick approach. [5l| 

• Viewing SOCs through the prism of conventional 
critical states has yet another flaw - it relies signifi- 
cantly on the renormalization group methods, which are 
essentially perturbative. On the other hand, perturba- 
tive methods can not be straightforwardly applied to 
avalanches (instantons), which are inherently of ultimate 
non-perturbative nature. The perturbation theory, how- 
ever, can be meaningfully applied to some collective vari- 
ables, (s^ - ls^ l In our case, such collective variables are 
the instanton modulii or rather the corresponding ghosts. 
The low-energy effective theory for the ghosts would pro- 
vide instantons with a sort of dual fermionic description, 
which admits perturbative treatment. This may be an 
advantage of the W-TFT picture of SOCs. 

• Soon after the introduction of the concept of SOC, it 
was proposed that SOCs may occur only in the so-called 
conservative models. [3, [Hal Whether the conservative- 
nonconservative classification is accurate is still under de- 
bate. [11113 The proposed W-TFT picture of SOCs did 
not rely on a specific form of the SDE. Hence, it may turn 
out that the SOC family is bigger that it is believed now. 
Some order-out-of-chaos-type systems [11] {e.g., pattern 
formations) may as well belong to the SOC family. 

• The origin of the self-tuning property of the SOC 
rests on the necessity to use the topological completion 
of the stochastic quantization (Parisi-Sourlas-Wu quanti- 
zation) due to the slow noise. This can be paraphrased in 
an informal yet seemingly accurate form: the criticality 
of SOC is "topologically" protected. 

• An important ingredient of all the W-TFTs is a set of 
Q-closed operators, from which one can construct topo- 
logical invariants ([2^ . Typically, these operators are cy- 



cles of different dimensionality in the base manifold and 
their explicit form depends on the field content of a the- 
ory and/or on topology of the target space. 

The fact that the Q-symmetry is broken may lead to 
the conclusion that these operators are not topological in- 
variants anymore and thus are useless. This is not quite 
so. The point is that due to the localization principle, 
perturbative vacua are still Q-closed (at least up to one- 
loop). Therefore, one can use these operators for the 
topological classification of the perturbative vacua. Hav- 
ing constructed such a classification, one may introduce a 
refined SOC time-evolution of the system as a sequence of 
jumps, at which the topology of the quantum state of the 
system suddenly changes. In general, however, not every 
physical avalanche changes the topology of the state. 

• To our opinion, one of the most interesting directions 
of future investigations could as well be the search for a 
possibly existing connection between SOCs and the con- 
cept of fault-tolerant (topological) quantum computing. 
[53 ] If exists, this connection would be an exciting news 
not only for physicists but also for neuroscientists. 

• We find it very unorthodox, though adequate to think 
of earthquakes as of fermions. 



VII. CONCLUSION 

In conclusion, let us sketch again the discussed sce- 
nario, which we believe is a likely candidate for a the- 
ory of a generic SOC behavior, and explicitly reveal the 
connections between the conditions for and properties 
of SOCs outlined in the Introduction. The condition 
of the slow external driving (C2) leads to the neces- 
sity of using the Stratonovich interpretation of noise in 
the stochastic differential equation representing an SOC 
(CI). This necessitates the topological "completion" of 
the stochastic quantization procedure (Parisi-Sourlas-Wu 
procedure) and leads to a model with BRST-exact action 
- to a Witten-type topological field theory. In the SDE's 
parameter space, there are full-dimensional regions, in 
which the BRST-symmetry is spontaneously broken and 
the SOC dynamics is liberated. The low-energy, long- 
wavelength part of the liberated SOC dynamics repre- 
sents avalanches (instantons) (PI). The liberated SOC 
dynamics can also be viewed as the Goldstone ghosts, 
which have no gap and thus are responsible for the criti- 
cal avalanche distribution (P2). The BRST-symmetry of 
the model and its breakdown, which are the essence of 
the criticality, can not be lifted by a moderate variation 
of the parameters of the model. This is the essence of 
the self-tuning property (P3). 

To the best of our knowledge, so far Witten-type topo- 
logical field theories have only been of " internal" math- 
ematical use as a tool for the studies of topologies of 
lower-dimensional manifolds. The proposal of this pa- 
per suggests that Witten-type topological field theories 
may find their applications in many other areas of science 
that study various realizations of self-organized critical- 
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ity, e.g., in geophysics, astrophysics, neuroscience, evolu- 
tionary biology, etc. 
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